The Efimov effect for three interacting bosonic dipoles 
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Three oriented bosonic dipoles are treated using the hyperspherical adiabatic representation, 
providing numerical evidence that the Efimov effect persists near a two-dipole resonance and in a 
system where angular momentum is not conserved. Our results further show that the Efimov features 
in scattering observables become universal, with a known three-body parameter, i.e. the resonance 
energies depend only on the two-body physics, which also has implications for the universal spectrum 
of the four-dipole problem. Moreover, the Efimov states should be long-lived, which is favorable for 
their creation and manipulation in ultracold dipolar gases. Finally, deeply-bound two-dipole states 
are shown to be relatively stable against collisions with a third dipole, owing to the emergence of a 
repulsive interaction originating in the angular momentum nonconservation for this system. 
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The recent advances in producing ultracold dipolar 
gases have sparked a great deal of interest in the novel 
phases that could be experimentally accessible [1]. Two 
key ingredients are the long-range and anisotropic nature 
of the dipolar interaction, which can be manipulated by 
applying an external electric field. This reorients current 
research scenarios in ultracold gases and opens up the 
possibility of intriguing phenomena across a broad range 
of different fields, from condensed matter physics to ul- 
tracold chemistry. Recent experiments [2] have been able 
to create a gas of ground-state dipolar molecules and ob- 
serve the electric field and geometry dependence of chem- 
ical reactions [3]. Nevertheless, a number of interesting 
effects can be expected when the external field is tuned 
near a dipole-dipole resonance (i.e., zero-energy dipolar 
bound state); this causes the interaction to vary from 
strongly repulsive to strongly attractive, in close analogy 
to the control through magnetic Feshbach resonances [4] . 

For nonreactive ground-state dipolar molecules [5] , in- 
elastic three-dipole scattering will be the main loss mech- 
anism, and these will determine the stability of the dipo- 
lar gas. A recent paper by Ticknor and Rittenhouse [6] 
predicts in fact that three-body recombination will be ex- 
tremely important for bosonic dipolar gases. Three-body 
recombination in this context, D -\- D -\- D D2 -\- 
involves three free dipoles colliding to form a dipolar 
dimer releasing enough kinetic energy to make the colli- 
sion partners escape from typical traps. As Ref. [6] points 
out, in the strongly interacting regime universal aspects 
related to the Efimov physics tend to dictate the colli- 
sional aspects of such processes. Nevertheless, the per- 
sistence of the Efimov effect for long range, anisotropic, 
dipolar interactions remains a open question. The Efi- 
mov effect [7] has proven to have a profound impact in 
both nuclear and atomic matter [8, 9]. In systems with 
short-range isotropic interactions, the Efimov effect man- 
ifests when the two-body scattering length as is much 
greater than the range of the interactions ro, through a 
series of three-body bound states with energies given by 
En^i/En = 6"^"^/^° (n = 0, 1, 2, ...), where the universal 



constant so= 1.00624 for identical bosons. Note that low- 
lying energies, however, are not universal and typically 
depend on the details of the interactions. In ultracold 
gases, the Efimov effect is explored near a Feshbach res- 
onance through its impact on scattering processes [10]. 

In this Letter, we show that the Efimov effect persists 
for dipolar systems, and moreover, that the dipolar in- 
teraction is extremely beneficial for the study of Efimov 
states. In particular, we find that the positions of the 
Efimov resonances (e.g. in energy, in ag, etc.), as well as 
other scattering properties, are universally determined by 
the two-dipole physics alone. Moreover, as the strength 
of the dipolar interaction increases, Efimov states tend 
to be increasingly long-lived. This scenario, introduced 
by the dipolar interaction, makes dipolar gases ideal sys- 
tems for the creation and manipulation of Efimov states. 
Following the method introduced in Ref. [11], we have 
also derived the scaling laws for three-body recombina- 
tion. Finally, we show the existence of an effective repul- 
sive long-range interaction between a dipolar dimer and 
a free dipole, which can prevent the decay of such dimers. 

The hyperspherical adiabatic representation [12] trans- 
forms the 3-body Schrodinger equation into coupled hy- 
perradial equations, given (in a.u.) by 
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where the hyperradius R describes the overall size of the 
system, and u represents all quantum numbers that la- 
bel each channel. In the above equation, fi = m/\/3 is 
the three-body reduced mass, E is the total energy, and 
Fi, is the hyperradial wave function in channel u. Pj^j^' 
and Quu' are nonadiabatic couplings that drive inelastic 
transitions, while the Uj^{R) are effective potentials that 
support bound and quasi-bound states of the system. 

The most stringent challenge in this approach is the nu- 
merical solution of the hyperangular adiabatic eigenvalue 
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FIG. 1: (color online) (a) A typical set of adiabatic hyperspherical potentials Uu{R) for three dipoles that exhibit the Efimov 
hyperradial potential curve (d^/ro = 38.11 and \as\ oo). Inset: rescaled adiabatic potentials showing the Efimov potential 
(the dashed lines are illustrations of the Efimov states in that potential, not to scale), (b) Schematic representation for the di 
dependence of as (upper part) and the Efimov spectrum (lower part) for dipolar systems (see text), (c) The Efimov potential 
signature at the first pole of as [di/ro — 4.86). The horizontal solid line shows the Efimov behavior given by Eq. (5) (see text). 



equation at fixed values of R in order to determine the 
potentials Uu{R) and channel eigenfunctions ^i,{R;Q)^ 
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^,(i?;l]) = 0. (2) 



Here Vt represents the set of all hyperangles describing 
the system's internal motion, is the grand angular 
momentum operator [12] and V = '^(ri2) + '^(r3i) + '^(7^23) 
is the pairwise sum of two-body interactions, given by 

v{r) = VosecHr/rof + ^^—^lllf^r). (3) 

The first term above is an isotropic short-range contri- 
bution and the second term is the anisotropic dipole- 
dipole interaction [with a short-range cut-off /(r)] be- 
tween dipoles aligned in the z direction with dipole length 
defined as (i^ = m <i^/2, where dm is the electric dipole 
moment. We emphasize that although our model as- 
sumes point dipoles it is also applied to ground-state 
dipolar molecules, where all the details of its compli- 
cated structure are encapsulated in the short-range be- 
havior of the interactions. Moreover, since the dipoles are 
free to interact in any geometry, the solutions we obtain 
from Eq. (2) already include both attractive and repul- 
sive aspects of the dipolar anisotropic interaction. To 
test whether any property of interest is universal with 
respect to the three-body short-range physics, we have 
performed calculations for different values of di as well 
as for different isotropic short-range potentials. 

The major difficulty introduced by the dipolar inter- 
action is that the three-body total angular momentum J 
is not conserved. To calculate Ujj{R)^ we expand ^ as 

!^) = ^ <Pi''{R; ip)DiM{a, p, 7), (4) 

J,K 



where I^^M Wigner D functions, a, /3 and 7 are 

the Euler angles, and (p are the Smith- Whitten hyper- 
angles [13]. Here M is the space- fixed frame projection 
of the total angular momentum J, K is the quantum 
number for its body-fixed frame projection, and H is the 
total parity. Finally, the body frame components at 
fixed R are numerically solved by expanding the and cp 
dependences in a B-spline basis. 

Here, we study the three-dipole problem for = 0^, 
which includes the = 0+ symmetry where the Efi- 
mov effect occurs for nondipolar systems. Evidently, for 
d£ <C To the Efimov effect might be expected to remain 
unchanged. However, for d£ > tq the coupling among 
different J's might possibly compromise its persistence. 
In order to investigate the Efimov effect for dipoles, we 
tune di > ro in Eq. (3) to get a two-dipole zero-energy 
bound state, where \as\ ^ ro [14]. Figure 1 summarizes 
our findings and sketches the three-body energy spectrum 
obtained. Figure 1(a) shows the adiabatic potentials Uu 
for di/ro = 38.11 and \as\ 00, showing the overall 
topology of the potentials, including a family of channels 
converging to deeply-bound dimers and channels describ- 
ing the collision of three free dipoles. For each dimer 
channel there exist an infinity of channels converging to 
that same threshold due to the coupling of all different 
J's [we have truncated the J-expansion (4) at Jmax = 14]. 

The inset of Fig. 1(a), however, shows the evidence 
of our central result. It confirms the existence of the 
universal Efimov potential. 
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which support the usual infinity of three-body Efimov 
states, illustrated schematically in the figure as horizontal 
dashed fines. Our numerical explorations suggest that, 
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14.534 
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25.498 


3.03x10-2 


6.6x10"^ 


38.110 


2.95x10"^ 


3.2x10"^ 


al^ /di ^ -8.1 


al^ /di ^1.8 


ahd/di ~ 8.6 



TABLE I: The positions and the widths of the lowest Efimov 
state for different values of di and the universal ratios for 
the positions of the Efimov features in three-dipole scattering 
observables. 



in the strong dipole limit where ^ ro, the rescaled 
potential in the inset of Fig. 1(a) is universal, i.e., the 
rescaled potentials for different values of di or for differ- 
ent strengths of the isotropic interaction, coincide almost 
perfectly. Figure 1(b) shows schematically, the three- 
body energy spectrum versus dipole length. As di > ro 
increases, new dimers are created and so are new fami- 
lies of Efimov states. The binding energy of the lowest 
Efimov state scales in proportion to l/mdj, based on our 
adiabatic potential in the inset of Fig. 1(a) as well as Ta- 
ble L Figure 1(c) gives further evidence for the long range 
universal Efimov potential when \as\ oo. We checked 
that adiabatic corrections to the Efimov potential are 
negligible. 

The existence of the Efimov effect for three dipoles 
might have been expected as a consequence of the 5- wave 
dominance of the two-body physics near the pole of ag [6]. 
Nevertheless, the Efimov effect for dipoles has a striking 
difference from the usual Efimov effect. As shown in the 
inset of Fig. 1(a) [see also Fig. 1(c)], the Efimov poten- 
tial extends only from values of R > di while for R < di 
the dipolar interaction causes the Efimov channel to be 
repulsive. This critically affects the universal properties 
associated with the dipolar Efimov effect. Since the Efi- 
mov states are well separated from the short-range region 
(at R ^ ro) by the repulsive barrier at < d^, properties 
of Efimov states are expected to be fully universal, i.e., 
they will depend on the two-dipole physics alone, namely, 
as and di. This implies that, in contrast to the usual Efi- 
mov effect, the three-body phase (or parameter) [9, 11] 
is now universal and one can derive the energies of the 
Efimov states (resonances) for both ag > and ag < 
{\as\ > di). Moreover, the barrier for R < di also sup- 
presses decay of the Efimov states and, therefore, the 
Efimov states for large values of di are likely to be more 
long-lived than in the nondipolar scenarios. 

To quantify the three-body parameter for dipolar in- 
teractions, we calculate the position and the width of 
the lowest Efimov resonance when \as\ ^ oo. To han- 
dle the sharp crossings between the Efimov channel and 
other deeply bound channels [see Fig. 1(a) and inset], we 
solve Eq. (1) by the slow variable discretization (SVD) 
method [15], and extract the position and width of the 



Efimov resonances. In Table I, we list the ground Efi- 
mov state energies Eq and the widths T for a few values 
of di. As we have expected, the positions of the Efi- 
mov resonances in Table I exhibit a universal trend. For 
the widths of these resonances, however, we do not ob- 
serve nor expect a purely universal behavior, since they 
encapsulate the decay to deeply-bound states at smaller 
distances. Nevertheless, for increasing di, the resonance 
width is suppressed as 1 /d^ , indicating increased lifetimes 
of the Efimov states. The universal three-body parame- 
ter n^=^/mEQ [9] is then estimated to be O.lT/d^ in the 
limit of di ^ ro- This knowledge of [9] then allows 
us to determine the universal formulas that predict some 
important three-body scattering observables. 
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Here is the rate coefficient for three-body recombina- 
tion into weakly (a^ > 0)- and deeply {as < 0)- bound 
dipole-dimers. a^d is the dipole plus dipole dimer scat- 
tering length, and Kei is the corresponding relaxation 
rate for such inelastic collision processes, D2 -\- D ^ 
D2-\-D. T] is related to the probability of decay to deeply 
bound channels [9]. Equations (6)- (9) now predict the 
positions of the main features of such rates that char- 
acterize the Efimov physics. For instance, the minima 
in Ks [Eq. (6)] should occur for values of given by 
4+ = l.Se^^/^od^ (n = 0,1,2,...). The ratio a^'/di 
determines the position of the Efimov resonances for Ks 
(as < 0) [Eq. (7)] while a}^^/di determines the position of 
dipole plus dipole-dimer Efimov resonances [Eqs. (8) and 
(9)]. The second part of Table I lists these properties in 
terms of ratios for all processes (we have dropped the fac- 
tor e^^/^° for simplicity). The universality on the three- 
dipole problem also has key implications for its four-body 
analog. In particular, the positions of the four-dipole fea- 
tures [16] can be obtained from our results in Table I. 

Our numerical results show that the attractive R~^ 
Efimov potential is absent when \as \ ^ di. Nevertheless, 
the adiabatic potentials are still universal. As shown in 
Fig. 2(a), the diabatic potentials show universal scaling 
with di. Figure 2(a) shows the lowest three-body con- 
tinuum channel and the lowest dipole plus dipole-dimer 
channel associated with the most weakly-bound dimer 
whose binding energy is proportional to 1/mdj. For 
this case, recombination to such dimer states proceeds 
through an inelastic transition near R ^ di, leading to 
the Ks oc dj scaling for recombination, consistent with 
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FIG. 2: (color online) (a) The diabatized hyperspherical po- 
tentials for the most weakly-bound dipole plus dipole-dimer 
channel with different values of di and \as\ < vq. (b) The di 
dependence of the asymptotic angular momentum lod and the 
expectation value of the angular momentum (p) for a deeply- 
bound dipolar dimer and the total angular momentum (J^). 



Ref. [6] . Since the dipole plus dipole-dimer channel is re- 
pulsive for R < di^we expect that recombination to more 
deeply-bound states will be substantially less important 
since the pathway to this final state would require tun- 
nehng through to small R < d^. 

The behavior of the dipole plus dipole-dimer channels 
when the dimer is deeply-bound also shows results in- 
teresting for dipole and dipole-dimer mixtures. These 
channels should be asymptotically described by 



W,{R) Ed- 



(10) 



where Ed is the dimer energy and lod represents its av- 
erage or effective orbital angular momentum relative to 
the free particle. For three particles with isotropic short- 
range interactions, l^d obeys the usual rules of addition 
of angular momentum. For instance, if the two-body 
state is an s-wave state, lod = J- For dipolar states, 
however, our results reveal that is no longer true and 
we rationalize it by the fact that the two-body state 
is not a pure angular momentum eigenstate. In fact, 
our numerical calculations show that lod becomes (in- 
dependent. Figure 2(b) shows that lod increases with d£ 
in the same trend as the expectation value of the angular 
momentum (j^) for the dipolar dimer and the expecta- 
tion value of the total (three-body) angular momentum 
(J^). However, the d^-dependence of lod is not universal, 
as it varies with the short-range physics. Nevertheless, 
the non-zero l^d provides an average centrifugal barrier 
that suppresses the dipole plus dipole-dimer relaxation, 
Vrei oc k'^Dd^r^''^^^ [11], where kod < ^o"^ the relative 
wave number. This suppression opens up the possibil- 



ity of creating stable mixtures of dipoles and two-dipole 
molecules. Further, with the tunability of the barrier 
that can be achieved via field control of d^, the interac- 
tion between a dipole and a two-dipole molecule can be 
controlled, providing an innovation for studies in ultra- 
cold quantum gases at both the few-body and many-body 
level. 

In summary, we have characterized the Efimov effect 
for three bosonic dipoles. A key conclusion is that long- 
range, anisotropic dipolar interactions are predicted to 
make the Efimov physics more universal than in the tra- 
ditional systems where this effect has been studied and 
observed. Efimov resonances with dipolar interactions 
are also predicted to be long-lived, which should aid their 
experimental observation and manipulation. The dipo- 
lar interactions introduce a tunable effective repulsion 
between a dipole and a dipolar dimer, which can be used 
to control the collisions between the molecules and also 
to stabilize the ultracold quantum gas of dipolar dimers. 
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